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Abstract 

We study how the synaptic connections in a pair of excitable electronic neurons affect the coherence of their spike trains 
when the neurons are submitted to noise from independent sources. The coupling is provided by electronic circuits which 
mimic the dynamics of chemical AMPA synapses. In particular, we show that increasing the strength of an unidirectional 
synapse leads to a decrease of coherence in the post-synaptic neuron. More interestingly, we show that the decrease of 
coherence can be reverted if we add a synapse of sufficient strength in the reverse direction. Synaptic symmetry plays an 
important role in this process and, under the right choice of parameters, increases the network coherence beyond the value 
achieved at the resonance due to noise alone in uncoupled neurons. We also show that synapses with a longer time scale 
sharpen the dependency of the coherence on the synaptic symmetry. The results were reproduced by numerical 
simulations of a pair of synaptically coupled FitzHugh-Nagumo models. 
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Introduction 

Neurons are highly nonlinear dynamical systems which are 
typically connected to tens of thousands of other neurons, the 
whole system being subjected to fluctuations whose stochasticity 
cannot be dismissed. This interplay between nonlinearity, high 
dimensionality and noise is what renders the brain a difficult and 
interesting system to study [1,2]. More generally, the last decades 
witnessed a surge in theoretical studies of collective phenomena of 
interacting nonlinear units. Since the seminal work of Kuramoto 
[3], for example, several aspects of synchronization have been 
addressed [4] . With the emergence of complex networks becoming 
a research topic in itself [5], the effects of topology on 
synchronization have been thoroughly investigated (see e.g. [6— 
8] for recent examples, or [9] for a review). Recendy, even the 
notion of networks of networks have emerged in the context of 
climate studies [10]. 

Even single neurons, however, can reveal surprises. In 1997, for 
instance, Pikovsky and Kurths unveiled the phenomenon of 
coherence resonance (CR), whereby an excitable system driven by 
white noise produces a spike train whose regularity (or coherence) 
attains a maximum at some finite value of the noise intensity [1 1]. 
In the low-noise regime, the spike train approaches a Poissonian 
incoherent behavior with small firing rate, whereas in the high- 
noise regime incoherence coexists with a large firing rate. At the 
resonance, the spike train looks almost periodic, despite the fact 
that the system is in an excitable regime, not tonic. 

The collective effects of coupling on CR were subsequendy 
investigated. It was shown that global coupling, either by square 
pulses [12] or via electrical synapses (gap junctions) [13], can lead 



to network synchronization with strong coherence. It was also 
shown that a network of excitable elements can exhibit system size 
CR, where increasing the number of elements in system leads, at 
first, to an increase in global coherence, while very large networks 
have reduced global coherence [13,14]. In addition, chemical 
synapses were shown to be better at increasing global coherence 
than gap junctions, even when the analyzed network contained 
only two neurons [14]. This highlights the importance of the 
characteristic times introduced by the chemical coupling in the 
post-synaptic response. 

Our aim in this contribution is twofold. First and foremost, we 
depart from previous studies on CR in networks of model neurons 
in that all of them have focussed on the effects of some coupling 
among the units which was uniform across the network: Wang et al. 
have employed the same uniform coupling intensity among all 
pairs of Hodgkin-Huxley neurons [12]; Toral and Mirasso had a 
single variable to parameterize the gap junction conductance 
among their FitzHugh-Nagumo elements [13]; and Balenzuela 
and Garca-Ojalvo simulated Morris-Lecar systems with chemical 
and electrical synapses: in either case, the same coupling intensity 
that connected neuron i to neuron j was also applied from j to i. In 
other words, previous analyses of collective effects have dramat- 
ically reduced the dimensionality of coupling space by focussing on 
the single-axis projection of uniform coupling. This is a convenient 
choice to render parameter space scanning feasible, but not very 
realistic. In nature, mutually connected neurons most probably are 
not symmetrically coupled. Here we set forth to investigate what 
happens to the coherence of spike trains when the coupling among 
the units is not necessarily uniform. Specifically, we fully explore 
the space of coupling intensities by focussing on the analysis of the 
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simplest network, namely, a pair of neurons. We will show that: 1) 
synaptic asymmetry can deteriorate the coherence of a pair of 
neurons that were previously uncoupled; 2) synaptic symmetry 
leads the system to a more coherent state, compared with the 
uncoupled case; 3) the effects of synaptic symmetry on the system 
coherence is strongly dependent on the characterstic time scale of 
the synapses. 

Second, we go beyond numerical simulations by employing 
type-II-excitable [15] electronic neurons which are connected via 
electronic circuits that mimic chemical synapses (see Materials and 
Methods). These electronic neurons are set at the edge of their 
Hopf bifurcation, as previously described [16]. The use of such 
electronic circuits gives our results an additional degree of 
robustness due to the variability of the electronic components 
and influence of external effects (such as thermal noise), all of 
which contribute to a more biologic-like scenario of noise and 
heterogeneities. Finally, we briefly discuss the potential applica- 
tions of our results beyond neuroscience, in other experimental 
setups where CR has been observed. 



dynamic variables, V = 1.5F_ — 0.67K OM( . In the actual electronic 
neuron circuit this is achieved with the use of an analog subtractor 
[16,17]. 

When set near its Hopf bifurcation, the electronic neuron is 
excitable, and its noise-induced spike train can be described 
approximately by the first two moments of the inter-spike interval 
(ISI) distribution P(t p ) [11]. The incoherence of the spike train is 

• , , , n \/<tl>-<tp> 2 

characterized by the parameter R p = — — ^jy and attains a 

minimum as a function of the noise intensity, which is controlled 
by the resistance R„oise (Fig. IB). Further details about the 
electronic neuron circuit, its model and the noise generator can be 
found in [16]. 

In order to connect two such electronic neurons, we have 
employed the electronic synapse shown in Fig. ID. Due to the 
diode at its input, it is activated whenever the output V pre of the 
pre-synaptic neuron becomes positive {V^}, which happens only 
during a spike [16]. Applying Kirchoffs law to the circuit, one 
obtains the dynamics for the voltage V c at the capacitor: 



Materials and Methods 

Electronic neurons and electronic synapses 

Previously, we have introduced a FitzHugh-Nagumo-like 
electronic circuit which models neurons with type-II excitability, 
operating in the ms time scale and receiving a noisy input with 
controllable intensity [16]. The circuit diagrams for both the 
electronic neuron and the noise generator are shown in Fig. 1A 
and Fig. IB respectively. The equations of motion that model the 
behavior of the electronic neuron can be readily obtained from 
Kirchoffs law and a simple model for the dynamics of the 
operational amplifier [16]: 

^ = S sr sign [ V h - V out + (V a - V h )®(aV out - F_ )] , (la) 



dt 



R 3 C 



[PV 0Ul + yV in -V_] 



(lb) 



where © is the Heaviside function, V a and Fj are the operational 
amplifier supply voltages and a = R\/(R\ + R 2 ), fi = Rn/(Rn + R5) 
and y = i? 5 /(i? 4 + R 5 ) (see Fig. 1A). The time scale for the 
dynamical variable V M is controlled by the slew rate S sr of the 
operational amplifier (Fig. 1A), with a typical value S sr = 20V/ /is. 
The characteristic time scale 1 / (R3 C) controls the dynamical 
variable V- and is set to yield V mt as the fast variable and V- as the 
slow one. The input variable V{ n receives the sum of a constant DC 
signal Voc, a noisy signal V„ 0 i !e provided by the noise generator 
(independently for each electronic neuron) and the synaptic input 
gVc (see details below) which couples the two electronic neurons: 



V in =V DC +V„ oise +gVc 



(2) 



The phase space of an uncoupled neuron (i.e. for Eqs. la and lb 
with g = 0 in Eq. 2) is shown in Fig. 1C. The nullclines V out = 0 
and V- = 0 resemble those of the FitzHugh-Nagumo model (see 
below) and the single fixed point can have its stability changed via 
a Hopf bifurcation (controlled by the Vdc) which generates a limit 
cycle and puts the system in a tonic regime. The membrane 
potential undergoing a spike (like the one showed in Fig. 2A, for 
instance), is obtained through a weighted subtraction of both 



d ^= 1 [V a3 (t)-V c ] 
at x 

RaRgC 



(3) 



where the characteristic time is x = R ^ and the asymptotic 

value is V a0 (t)= + r^ J V p + re {t). The electronic analog V s of the 

excitatory post-synaptic potential (EPSP) is obtained by multiply- 
ing V c by a controllable gain with a standard amplifier [17] which 
effectively controls the strength of the coupling between the two 
electronic neurons: 



V s =gV c 



(4) 



Here we focus on excitatory synapses (g > 0), but an inhibitory 
synapse (i.e. with g<0) can easily be mimicked with an inverter 
amplifier. First-order kinetics such as that of Eq. 3 is considered a 
reasonable approximation for the dynamics of some classes of 
synapses (e.g. based on AMPA or GABA receptors) [2]. 

As shown in Fig. 2, this setup allows us to control several 
interesting features of the electronic EPSP: R x controls the rise 
time of V s (Fig. 2A). Both R^ and Rp have influence on the decay 
time of V s , as well as on its maximal value S (Fig. 2B), which can 
also be independendy controlled by the gain g (Fig. 2C). 

With a neuron in its excitable regime, we can control its 
spontaneous firing rate by adjusting the root mean square (rms) 
Vrms of the zero-mean noise voltage V„ 0 Lse that stimulates it (i.e. 

Vrms — \J *\ Vnoise))- Figure 2D shows a time series of a pre-synaptic 
neuron (labeled 1) with noise intensity ~ 1.2(1 )V and the 
corresponding EPSPs it generates whenever it spikes. A post- 
synaptic neuron (labeled 2), also in its excitable regime, receives 
these EPSPs as well as noise, but with lower intensity 
Vfms — 1-0(1)V (noise sources in different neurons are indepen- 
dent). Driven by this lower- intensity noise alone, the post-synaptic 
neuron is extremely unlikely to fire, which it eventually does upon 
receiving an EPSP with peak voltage Sn = 1.3 V (see Figs. 2D and 
E). In Fig. 2F we show that the post-synaptic spike probability 
increases monotonically with the synaptic amplitude S, with a 

sensitivity threshold that decreases with increasing noise intensity 
ypost 
rms ' 

In what follows, experimental results were obtained with 
electronic neurons whose parameters were chosen to be as similar 
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Figure 1. Electronic circuits. (A) Excitable electronic neuron circuit. V a = — Vt = 12.0V are the supply voltages, V oul and K_ are the dynamic 
variables and K,„ is an input voltage. Voc= — 7.30(2)V. (B) Noise generator circuit. The noise intensity is controlled by the resistance R, u ,i Se . (C) Phase 
space diagram for the system described by Eqs. 1 . The nullclines V mll = 0 (black solid line) and F_ = 0 (black dashed line) intersect at the stable fixed 
point. Due to the noise provided by the noise generator at the input V m , the system is forced out of the rest state, often having to perform a long 
excursion in the phase space, producing a spike (red line). (D) Circuit that mimics the behavior of a chemical synapse. V pn , is the pre-synaptic input, 
while V+ e takes into account only positive values of V pre . The resistances R x , Rp and the capacitor C= 10 nF control the time scale z of the synapse 
(see text for details). The equivalent of the synaptic conductance g can be set by a standard amplifier, which yields the post-synaptic potential 
V s =gVc (see text for details). 
doi:1 0.1 371 /journal.pone.0082051 .g001 



as possible (within the ~5% tolerance of the electronic compo- 
nents). Where not shown, error bars are smaller than symbol sizes 
and uncertainties in experimental values are in the range of 1%. 
The noise intensity will be denoted by the resistance R no ise 
controlling the gain of the noise amplifier, as the voltage rms 
increases linearly with R no ise [16]. Experimentally, the peak value 
(amplitude) S of the EPSP Vs (see Fig. 2B) is easier to measure 
than the gain parameter g, and will therefore be used as a measure 
of the synaptic strength, with Sy denoting the coupling from 
neuron i to neuron j (with (ij)e{(l,2),(2,Y)}, in the present case). 

Computer simulations 

Numerical simulations using the FitzHugh-Nagumo model for 
neuronal excitability were made to confront the experimental 
results. This model was chosen due to its simplicity and close 
similarity with our model for the electronic neuron dynamics (see 
Eqs. 1 and 2): 



dx 
' dt ' 



x 

T 



dv 

-^=x + C+D£(t)+gVc , 



(5a) 



(5b) 



where £ = 0.001 is set to reproduce the ratio of the time scales 
governing the variables in Eqs. 1 and C is an assumed delta- 
correlated Gaussian noise with zero mean whose intensity is 
controlled by the parameter D. Similarly to Eqs. 1 and 2, the 
strength of the synaptic coupling between the two model neurons 



is controlled by the gain constant g. The parameter £ = —1.05 
ensures that each FitzHugh-Nagumo model is in a excitable 
regime but very close to its Hopf bifurcation [11]. Coupling 
between the two FitzHugh-Nagumo systems is achieved using the 
model for the electronic synapse, as previously described in Eq. 3: 



~ = ~[0T 0(x(O) x(t)-v c 
at t 



(6) 



where again v c integrates pre-synaptic activity x only when it is 
positive. We will employ the gain g y as a measure of the coupling 
strength from model neuron i to model neuron j. The equations 
were integrated using Euler-Maruyama's method with a time step 
A? = 0.005. 

Results and Discussion 

Coherence deteriorates with an incoming synapse, but is 
restored with an outgoing synapse 

To investigate the effects of symmetry in the synaptic coupling 
on the coherence of spike trains, we started by the asymmetrical 
extreme of connecting two neurons unidirectionally. As shown in 
Fig. 3 A, the incoherence R} p of the pre-synaptic neuron 1 exhibits a 
minimum as a function of the noise intensity, as is typical of CR 
[11,16]. Due to the synapse from neuron 1 to neuron 2 (with 
synaptic amplitude S12 = 1-0V), the spike trains of neuron 2 are 
less coherent than those of neuron 1 (Fig. 3A). This could be 
expected, since neuron 2 is receiving noise-induced spikes from 
neuron 1 in addition to its own (independent) noise source. In this 
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Figure 2. Output V s of the electronic synapse when subjected to a pre-synaptic spike. (A) Black line corresponds to the pre-synaptic spike 
(scale on the right). The rise time of V s is governed by R„_ (Rp = lOkO, R g = 2kQ), while (B) changes in both R„ and Rp influence the decay time of V s 
(when V+ e (t) = 0V; R„ = 5kQ, R g =2kQ). (C) The gain g (see text for details and Fig. 1) is proportional to the amplifier resistance R s , which affects 
only the maximum amplitude 5 of Vs {R tl = 5kQ, Rp = 20kQ). (D) Spike train of a pre-synaptic neuron excited by noise (Vg£ 1.2(1)V) and the 
resulting Vs with S\2 =; 1.3V (red). (E) Spike train of a post-synaptic neuron (blue) excited by Vs and a lower noise intensity (K^f ~ 1.0(1)V). (F) Post- 
synaptic spike probability as a function of S for two values of K/JJJf. In both cases, the pre-synaptic neuron is subjected to V?* ss 1.2(1)V. 
doi:10.1371/journal.pone.0082051.g002 



simple scenario, the behavior of the pre-synaptic neuron is, as 
expected, unaffected by its outgoing synapse. 

Keeping the synapse from neuron 1 to neuron 2, we added 
another synapse in the opposite direction, with the same synaptic 
amplitude (S21 = 1 .0V). This reverted the effect of the single 
synapse, causing neuron 2 to reattain a minimum in the R p vs 
Rnoise curve, as shown in Fig. 3B. This is perhaps counterintuitive, 
since one might expect that, by synaptically coupling the less 
coherent neuron-2 spike train with neuron 1, R l p should increase. 
What happens, however, is that not only the coherence of neuron 
1 is weakly affected, but also neuron 2 recovers its coherence. 
More importantly, it does so by means of an outgoing synapse. 

We explored how the above phenomenon unfolds as we 
gradually change the synaptic strengths. We started with initially 
uncoupled electronic neurons (Si2 = S2i =0V) and noise intensity 
near its resonance value (R m i se = l0k£2). This choice of noise 
intensity was made in order to maximize the variation of the 
incoherence parameter R p as the coupling is varied (see Fig. 3A 
and Fig. 3B). Increasing only the synaptic strength S12, 
incoherence R 2 p in the post-synaptic neuron increased monoton- 
ically, while R p remained essentially unchanged (Fig. 3C). With 
£12 = 1.4V fixed, we then increased S21, which led to a rapid 
increase in the coherence of neuron 2. Neuron 1, on its turn, 



showed a small decrease of incoherence (Fig. 3D), in a 
phenomenon similar to what has been reported in numerical 
simulations of symmetrically coupled neurons [14]. 

Coherence depends weakly on synaptic symmetry for 
fast synapses 

The above results suggest that symmetry between the synaptic 
strengths Sn and S21 plays an important role in the spike train 
coherence of both neurons. To perform a thorough investigation 
of this phenomenon, we looked into the dependency of both R p 
and R p on S12 and S21 in a large region of the parameter space. In 
Fig. 4A we show the EPSP for both synapses with different 
synaptic strengths and the corresponding spike times on both 
neurons when the synaptic time scale is t = 0.04ms (in what 
follows, Ra = 5kQ is kept fixed and 1 is controlled only by the 
resistance Rp). The dependency of the incoherence parameter R p 
on the synaptic strengths is shown in Fig. 4B (of which Figs. 3C 
and 3D are cross-sections). The firing rate of neuron 2 was also 
measured as a function of the synaptic strengths and is shown on 
Fig. 4C. Note that in this case an increase in the excitatory synapse 
(with strength S21) from neuron 2 induces a decrease in its own 
firing rate, which leads to an increase in the spike train coherence. 
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Figure 3. Incoherence parameter R p as a function of the noise intensity and synaptic amplitudes. (A) Neuron 1 (pre-synaptic) is 
connected to neuron 2 (post-synaptic) through a synaptic circuit with parameters R% = 5kO, Rp = 20kO_ and synaptic amplitude Sn = IV. IJj presents 
a minimum, as expected, but neuron 2 has increased incoherence Rj r (B) Addition of an identical reciprocal connection leads to a coherence recovery 
in neuron 2. (C) R p and Rj as functions of the synaptic amplitude S u when S 2 \ =0V (unidirectional connection). (D) R p and R 2 p as functions of the 
synaptic amplitude S21 with 5i2= 1.4V. In (C) and (D), noise is set at the coherence resonance value for a single neuron (R„ 0 j sl , = 10kO). 
doi:1 0.1 371 /journal.pone.0082051 .g003 



The parameters R p and F\ for neuron 1 behave in the same way if 
the indexes 1 and 2 are exchanged. 

The total effect of the synaptic coupling on the system can be 
measured if we define the mean incoherence parameter 
(R p y = (R 1 p + Rp)/2. We show its dependence on Su and S21 in 
Fig. 4D. This graph suggests that the incoherence parameter can 
be described as a function of symmetry and asymmetry parameters 
defined as s = Su + S21 and a = Su — S21 respectively. The s and a 
axes are also shown in Fig. 4D. Although for large values of \a\ an 
increase in (R p } is observed, there is a large flat region in which 
no significant change in coherence is observed. The larger the 
symmetry s, the larger the flat region along the a axis (Fig. 4E). For 
fixed asymmetry a, increasing s causes an increase in the system 
coherence, as shown Figure 4F. 

Coherence depends strongly on synaptic symmetry for 
slow synapses 

The flat region in Fig. 4D might suggest that synaptic symmetry 
is not so relevant for coherence, after all. This scenario changes 
significandy, however, when the synaptic time scale increases. 
Figure 5 shows similar results to those of Fig. 4, but with 
t~ 0.048ms (ity = 100kn). This 20% increase in the value of the 
synaptic time, as compared with the previous value, leads to a 
significant qualitative change in the EPSPs, as shown in Fig. 5A. 

For these slower synapses, R 2 p now displays a narrow valley of 
local minima as a function of S12 and i?2l (Fig. 5B). If one revisits 
the situation in which Su # 0V and S21 is increased from zero, 
once more an increase in coherence is observed in neuron 2 owing 
to an outgoing synapse (Fig. 5B). Differently from the scenario of 
the fast synapses (Figs. 4B and 4C), however, now the increase in 
the coherence of neuron 2 occurs with an increase in its firing rate, 
as shown in Fig. 5C. 



The effect of synaptic symmetry on the overall coherence is 
much more pronounced for slower synapses, as shown in Fig. 5D. 
A much sharper minimum of (Rp} emerges near a = 0, regardless 
of the value of s (Fig. 5E). Furthermore, for fixed a — 0, an increase 
in symmetry j (which amounts to an increase in overall synaptic 
strength) leads to an increase in coherence, which attains values 
above those seen for uncoupled neurons (Fig. 5F). 

Comparison with computer simulations 

We attempted to reproduce the above results in numerical 
simulations using the widely known FitzHugh-Nagumo model, as 
described in Eqs. 5. As in the experiments made with the 
electronic circuits, we set the model neurons in the excitable 
regime ((=—1.05 in Eq. 0) but close to its Hofp bifurcation (at 
£= — 1), so that the system can fire under the influence of the 
Gaussian noise input £ with intensity D. The two model neurons 
are connected through our model of the electronic synapse (Eq. 6) 
with a coupling coefficent g (see Eq. 5). The resulting EPSP gvc 
generated by a pre-synaptic spike is then added (alongside with the 
noise) to the variable £ of the post-synaptic model neuron (see Eq. 
5b). We find the noise intensity that yields the minimum of the 
incoherence parameter R„ for the uncoupled case (Z) = 3 x 10~ 3 ) 
and then vary the synaptic strengths gu (from the model neuron 1 
to model neuron 2) and £21 (from 2 to 1) measuring the mean 
incoherence parameter <-Rp> for each value of the pair (g\2,gi\) 
as we did with the electronic circuits. Note that the approach here 
is to measure the synaptic strength direcdy from the amplification 
factors g\2 and £21 instead of measuring the amplitude of the 
EPSPs. The simulations were performed for two different values of 
the synaptic time scale T, first for T = 0.1 (fast synapse) and then for 
T= 1.0 (slow synapse). 

The results of the numerical simulations are shown in Fig. 6. 
Comparing first Fig. 6A with Fig. 4D, we observe that the 
numerical model corroborates the results of the electronic circuits: 
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Figure 4. Influence of symmetry on CR for fast synapses. Noise intensities are set at the coherence resonance value (R m ,i„ = lOkQ) for each 
neuron. Synaptic parameters are Rp = 20kf2 and x— 0.04ms. (A) Time series showing spike times for neuron 1 (black dots) and neuron 2 (blue dots) as 
well as the EPSP of the synapse connecting neuron 1 to neuron 2 with Si 2 = 0.4V and the reverse synapse with S%\ =0.2V. (B) Incoherence R 2 p of 
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horizontal plane. (E) {R p } vs the asymmetry parameter a. (F) {R,,} vs the symmetry parameter s. 
doi:10.1371/journal.pone.0082051.g004 



with fast synapses, the dependence of the incoherence on the 
synaptic symmetry is weak. The mean incoherence parameter as a 
function of the symmetry parameter j and asymmetry parameter a 
is shown in Fig. 6B and Fig. 6C and there is good agreement with 
Fig.4E and Fig. 4F. Note that in the case of the computer 
simulations, we employ U=gn — g2\ and S = gu+g2\- 

Next we increase the synaptic time scale of the model (Fig. 6D) 
and compare the results with those of Fig. 5D. Although the shape 
of the surface obtained through the simulations does not exactly 
reproduce the one obtained experimentally with the electronic 
circuits, the dependence of the mean incoherence parameter on a 
and s keeps its main features for the case of slow synapses: the 



(Rpy dependency on a is sharper than in the case of fast synapses 
(Fig. 6E). It also reaches smaller values when s increases, as 
compared to faster synapses (Fig. 6F). 

Conclusions 

We have studied the effects of synaptic coupling between 
excitable electronic neurons on the coherence of their spike trains. 
We have shown that the coherence resonance observed in 
uncoupled neurons deteriorates in the post-synaptic neuron when 
the synaptic connection is unidirectional. Furthermore, coherence 
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Figure 5. Influence of symmetry on CR for slow synapses. Same as Fig. 4, except that synaptic parameters are Rp = lOOkD and t~ 0.048ms). 
doi:1 0.1 371 /journal.pone.0082051 .g005 



resonance is restored and strengthened when the synaptic loop is 
closed in a bidirectional coupling. 

It is important to emphasize that, although this phenomenon 
may look similar to the synchronization between bidirectionally 
coupled spiking neurons [4, 1 8] , here neurons are excitable and the 
only attractor of the system is a fixed point with both neurons 
quiescent. In this sense, it is interesting that the interplay between 
noise and synaptic coupling leads a post-synaptic neuron to regain 
its coherence by means of an outgoing synapse. 

Furthermore, the mechanism by which this increase in 
coherence is attained depends on the synaptic time scale. With 
faster synapses, moving from a unidirectional coupling to 
bidirectional coupling by strengthening one of the synapses leads 
to an increase in coherence while firing rates decrease (despite the 
fact that all synapses are excitatory). For slower synapses, the same 
coherence increase is achieved with an increase in firing rates. 



Overall, our results show that, for fast synapses, the average 
coherence of the spike trains can be maintained in a broad region 
of synaptic-strength parameter space. However, as synapses 
become slower, maximal coherence is achieved only in a much 
more restricted region, around the symmetry axis Sn = S21 . Along 
this axis, strengthening synaptic connections lead to an increase in 
coherence beyond the values attained by isolated neurons. 

All the results above for electronic neurons were qualitatively 
reproduced with computer simulations of synaptically coupled 
FitzHugh-Nagumo models, suggesting the phenomenon is robust. 
It would be worth exploring whether it remains valid when 
neurons have a different excitability class, such as those near a 
saddle-node bifurcation [19]. Such type-I-excitable neurons can 
often be further reduced to simpler descriptions [20] which then 
might allow an analytical understanding of the results presented 
here. 
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0.1 0.15 0.2 



Figure 6. Influence of symmetry on CR for numerical simulations of synaptically coupled FitzHugh-Nagumo models. We have 
employed e = 0.001, J= — 1.05 and fl = 3x 10~ 3 in all simulations. For t = 0.1 (fast synapses): (A) Mean network incoherence (R p } vs g\% and g%\. The 
symmetry and asymmetry axes are drawn on the horizontal plane. (B) {R p } vs the asymmetry parameter a. (C) {R p } vs the symmetry parameter s. 
(D), (E) and (F) are the same as (A), (B) and (C), except that synapses are slower: t= 1.0. 
doi:1 0.1 371 /journal.pone.0082051 .g006 



Naturally, the ultimate test of our results would come from 
electrophysiological recordings of real neurons in which the 
symmetry of the synaptic coupling could be controlled. Although 
we are unaware of experiments in that direction, our results could 
also be useful to (or validated in) other systems not directly related 
to neuroscience, but where CR has been experimentally observed. 
The experimental setups range from semiconductor lasers [2 1] to 
the famous Belousov-Zabotinsky chemical reaction [22], recently 
also reaching nanoscopic scales in the ionic transport through 
single- walled carbon nano tubes [23]. If, for instance, applications 
on these setups depend on highly coherent states, then symmetry 
in the coupling may prove useful to overcome coherence 
degradation due to external factors. 

The emergence of self-sustained activity, a recurrent theme in 
the field of excitable media [24], could also be the subject of future 
studies regarding the effects of coupling asymmetry. Even in our 
simple two-neuron network, for example, preliminary results 
suggest that increasing the time constant % beyond the values used 
in this work can throw the system in an self-sustained attractor, 
despite the fact that both neurons are individually in a excitable 
state. Similar effects have been observed in a pair of electrically- 



coupled 0-neurons [20] and are likely to be important in an 
electronic implementation of neuro-inspired artificial sensors, 
which are predicted to have maximal dynamic range and 
sensitivity at the transition to a self-sustained state [25]. 

Taken together, our results point to the importance of allowing 
for inhomogeneity in CR studies of coupled excitable elements. 
Our study of a coupled pair, where inhomogeneity is reduced to 
the asymmetry of the synaptic connections, can be regarded as a 
first step towards larger networks. 
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